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Abstract 

In this paper we study the asymptotic behaviour of a nonlocal non- 
linear parabolic equation governed by a parameter. After giving the exis- 
tence of unique branch of solutions composed by stable solutions in sta- 
tionary case, we gives for the parabolic problem L°° estimates of solution 
based on using the Moser iterations and existence of global attractor. We 
finish our study by the issue of asymptotic behaviour in some cases when 
i — >■ oo. 
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1 Introduction 

The non-local issues are important in studying the behavior of certain physi- 
cal phenomena and population dynamics. A major difficulty in studying these 
problems often lie in the absence of well-known properties as maximum princi- 
ple, regularity and properties of Lyapunov (see [S], [B]) and also the difhculty 
to characterize and determine the stationary solutions associated thus making 
study the asymptotic behavior of these solutions very difficult. 
In this paper we study the solution u(t, x) to the nonlocal equation 

' Ut — div{a{lriu{t)))\7u) = f dans K+ x fl 

u{x,t)^0 sur M+ X 917 (1) 
0) = Mo dans Q.. 

In the above problem uq and / are such that 

u^^L\n), f €L\Q,T,L\n)), (2) 

with T a arbitrary positive number, a is a continuous function such that 

3m, M such that < m < a(e) < Af Ve e M. (3) 

The nonlocal functional Ir is defined such that 

lr{.){x):L'{n)^% u ^ lr{u{t)){x) = j g{y)u{t, y)dy . (4) 

nr\B(x,r) 

Here B{x,r) is the closed ball of K" with radius r and g G L^(r2).. It is some- 
times possible to consider g more generally, especially when one is interested in 
the study of stationary solutions of (see |3]). 

In physical point of view problem ([T]) gives many applications especially 
where g = 1 in population dynamics. Indeed, in this situation u may represent 



a population density and lr{u) the total mass of the subdomain 17 n B{x, r) of 
17. Hence ([1]) can describe the evolution of a population whose diffusion velocity 
depends on the total mass of a subdomain of 17. For more details of modeli- 
sation we refer the reader to [3- This type of equation can be applied more 
generally to other models including the study of propagation of mutant gene 
(see [H] , [H] , ) ■ A very recent study of this propagation was made by Ben- 
dahmane and Sepulveda [5] in which they analyze using a finite volume scheme 
adapted, the transmission of this gene through 3 types of people: susceptible, 
infected and recovered. 

In mathematical point of view, when r — d where d is the diameter of 17 
problem ([1]) has been studied in various forms (see f6l.[8].[9].|15|). 

However when < r < d, several questions from the theory of bifurca- 
tions have arisen concerning the structure of stationary solutions including the 
existence of a principle of comparison of different solutions depending on the 
parameter r and the existence of branches (local and global) of solutions. A 
large majority of these issues has been resolved in [3]. It shows that when a is 
decreasing the existence of a unique global branch of solutions and existence of 
branch of solutions that are purely local. Some questions may then arise: 

(i) The unique branch described in [3] it is composed of stable solutions? 

(ii) What about stability properties of the corresponding parabolic problem? 

The plan for this work is the following. In section 2 we give some existence and 
uniqueness results. Section 3 is devoted to stationary problem corresponding to 
((!]). In particular, we study in a radial generalisation of Chipot-Lovat 

results about determination of the number of solutions. We also establish that 
the unique global branch of solutions described in [3] is composed by stable 
solutions (theorem l3.8p . In section 4 firstly we address a L°° estimate taking to 
account estimate based on Moser iterations. Secondly we prove existence of 
absorbing set in , which allows us to prove the existence of a global attractor 
associated to ([T]) (see remarkjS]) . Finally we obtain a result of stability properties 
of the corresponding parabolic problem. 

2 Existence and uniqueness results 

In this section we show a result of existence . We set V = Hl{Vl) and V' its 
dual, we take the norm in V, \\.\\v such that 



< . , . > means the duality bracket of V' and V. 
Then we have 

Theorem 2.1. LetT>0,fe L'^{0,T,V') and uq G L^(17), we assume that a 
is a continuous function and the assumption checked then for every r fixed, 




n 



2 



r G [0, diami^Q,)], there exists a function u such that 
we 2.2(0, T,F), Ut e L^{0,T,V') 



it(0, .) — uq in il 



f^{u,(f>) + J^^a{lr{u{t)))Vu\7(pdx=< f,(f>> in D'{0,T) \/(P e H^{n). 



(5) 

Moreover if a is locally Lipschitz i. e 

Vc such that |a(e) — a(e')| < 7c|e — e'| Ve, e' G [— c, c], (6) 

then the solution of is unique. 

Remark 1. Before to do the proof it is necessary to see that for r — problem 
^ is linear and the proof follows a well-known result (see 110^). it is even 
when r = diam{n) (see We will focus therefore in the following where 

r g]0, diam{D,)[. 

Proof. For the existence proof we will use the Schauder fixed point theorem. 
Let w G i2(o,T,i2(f7)) we get 

t-^lriw{t)), 

is measurable as a is continuous then 

t—^a{lr{w{t))), 
is too. The problem of finding u — u{t, x) solution of 

u G L^{0, T, V) n C([0, T],L^{n)) Ut G £2(0, T, V) 
u{0, .) = uo 

f^{u,(j)) + ^^a{lr{w{t)))'sJuV(j)dx =< f,(j)> in D'{0,T) Vc/) e H^ifl), 

(7) 

is linear, besides ([7]) admits a unique solution u = Fr{w) (see[TU],[7]). Thus we 
show that the application 

w — >Friw)=u, (8) 

admits a fixed point. Taking w = u in ([7]) we get using ^ and the Cauchy- 
Schwarz inequality 

ll\u\l + m\\u\\l.<\fUMv, (9) 
is the usual norm in V and |/|^ is the dual norm of /. We take 

T 1 
2 



\u\l^q,t,v} = \ / \\u\\vdt 







Using Young's inequality to the right-hand side of it follows that 
By integrating (fTO|) on (0, t) for t < Twe obtain 



t 



Iwm + y / hfydt < l\uo\i + ml (11) 



3 



We deduce that there exists a constant C — C{m, uq, f) such that 

\u\mo,T,v) < C (12) 

Moreover 

< ut,v > + < -~div{a{lr{u(t)))Vu), v >=< f,v > Vv e V, 
This gives us 

\utU<M\\u\\v + \fU. (13) 
By raising (|13p squared and using the Young inequahty we have that 

\ut\l<2M^ufy + 2\f\l. (14) 

By integrating ([T^ on (0,t) and assuming (TT^ we obtain 

|wt|L2(o,T,y') < C*', (15) 
with C" = C"(m, M, /, ug) and C is independent to w. It follows from (fT2| and 

m 

l^tlL2(o,T,y') + I"Il2(o,t,v) ^ (16) 
with i? = + C'^ . From (fT2)) and the Poincare inequality it follows that 

\u\L^{o.T.L^{n)) < R', (17) 

By setting 

Ri=max{R',R), (18) 

and associating ([T7)) and ([T5)) . it follows that the application F maps the ball 
B{0, Ri) of L^{0, T, L^{n)) into itself. Moreover the balls of ^^^(O, T, V, V) are 
relatively compact in L^(0, T, L^(f2)) (see for more details), ([TB]) clearly 
shows us that F{B{0, Ri) is relatively compact in B{0,Ri) with 

B{0, Ri) = {ue L\0, T, L\n)); \u\L^o,T,L-m) < Ri}- 

In order to apply the Schauder fixed point theorem, as announced, we just need 
to show that F is continuous from B{0,Ri) to itself. This is actually the case 
and completes the proof of existence. 

We will now discuss the uniqueness assuming of course that assumption ^ 
be verified. Consider ui and U2 two solutions ([5]), by subtracting one obtains in 
D'{0,T) 

— (wi-U2,w)+ / {a{lriui{t))Vui{t)~a{lr{u2{t)))\/u2{t))\/(j)dx = y<f)eH^{n). 
n 

(19) 

Since 

a(Zr(wi(t)))VMi - a{lriu2it)))Vu2{t) = {a{lr{ui{t))) ~ a{lr{u2{t)))Vui{t) 

+ a{lr{u2{t)))V{ui{t) - U2{t)), 

(20) 
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we get 

d f 

— {ui - U2,v)+ / a{lr{u2{t)))V{ui{t) ~ U2{t))\7(pdx 
n 

^ - J {a{lr{ui{t))) - a{lr{u2{t)))VuiS/(j)dx V(/) £ H^{^). 
n 

(21) 

Moreover ui,U2 G C{[0,T], L'^{fl)) there exist z > such that 

lr{ui{t)),lr{u2{t)) e[~Z,z]. (22) 

Taking v — ui — U2 ui (I2ip . it comes easily by Cauchy-Schwartz inequahty and 

m 

- U2\l + m\\ui - U2\\v < 7l'r(wi(i)) - lriu2it))\\\ui\\v\\ui -U2\\v- (23) 
We get in [3] 

\lrHt)) < C\Bix,r) nn\^-\g\2\uit)\2 < \n\^\g\2\u{t)\2, (24) 

where C a constant, |r2| represents the measure of Q and n V 3 the maximum 
between the dimension n of and 3. By using p4L ((23)) and the Young in- 
equahty 

1 ,9 m , 

ab < H a . 

- 2m 2 

We deduce 

— \ui - U2I2 + "T-llui - U2\\y < Pit)\ui ~ U2I2, (25) 

with 

p{t) = 1(^C|1]|^|5|2 ll^ilk)" e ^'(0,T), 



which leads to 



^\ui - U2\l < p{t)\ui - U2\l- (26) 



Multiplying dM]) by e" ■/'o it follows that 

e-^'P^^^'^^^lui -U2I2 -p(t)e-foP^'^'^'\ui-U2\l < 0. (27) 

Hence 

l{e-/oP(^)'i-|ui-7i2|2}<0. (28) 

t 

— J p{s)ds 

This shows that t 1 — >■ e " |ui — U2I2 is nonincreasing. Since for t = 0, 

ui(0, .) = 1*2(0, .) = Uq. 

This function vanishes at and nonnegative, we conclude that it is identically 
zero. This concludes the proof. □ 
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3 Stationary solutions 

Consider the weak formulation to the stationary problem associated to ([T]) 
rp^i~div{a{lr{u))Vu)^f dans il 

3.1 The case r = d 

By taking cj) the weak solution of the problem 



-A0 = / dans 51 



(30) 



we get due to a Chipot-Lovat [8 results that 

Theorem 3.1. Let a be a mapping from K. into (0, oo). The problem (Pd) has 
many solutions as the problem in M 

lia{n) = ld{(j)), (31) 

with /I = Idiud)- 

Remark 2. Theorerr A3.1\ allows us to see where a is increasing that the problem 
Pd admits a unique solution and determine for a given a the exact number of 
solutions (Pd)- However it is difficult or impossible to adapt the proof of the 
theorerr iS.ll in case < r < d. 

3.2 The case < r < d 

As announced in the introduction we focus our study to the case of radial 
solutions of (Pd)- We will assume is the open ball of with radius d/2 
centered at zero. We set 

Ll{n) = {u e L^{n) 3u e L2(]0, d/2[) such that uix) ^ u(|la;||)}, 

and we also assume that 

feLliQ) 

aeW^'°°{R),Ma>0 (32) 

/ > a.e in H. 
g > a.e in fl. 

We start by giving in some sense in a linear case a result that will be used 
later to explain the asymptotic behavior. 

Proposition 3.2. Let A,Be C{fl) be positive radial functions such that A < B 
in n and also f,h E L'^{fl) two positive radial functions. Let u G H^{Q)) the 
radial solution to 

- div{A{x)Vu) = f in ft, (33) 

and 

- div{B{x)Vu) = h in Vt. (34) 

Then f < h a.e in fl. 
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Proof. We proved in [3] that if u is a the radial sohition of ([55)1 then for a.e t in 
[0,d/2], 

t 

u'it) = -J^Jj'-r~\f{s)ds. (35) 

From dMl) and ^ we obtain 

t t 

B{t) 



□ 





Since A < B in H and /, > with / ^ 0, /i ^ hence f < g. 



In a nonlocal case, some results of existence of radial solutions and compar- 
ison principle between Ur, Ud and uq has been demonstrated in [3]. It is also 
proved that if we set for all r G [0, d] 

/,. := [inf ('/'), sup ^r(0)]- (36) 

Here denotes the solution of 

f-A0==/ dans 17 

\0ei/i(r!). 

By the inclusion or not of Ir at an interval of R we somehow generalize the 
theorem l3.1l 

Lemma 3.3. Let r G [0,c?]. Assume that hS'A) holds true and there exist < 
mi < 7712 suc/i t/iai 

a(mi) = max a a{m2) = min a (38) 

[mi,m2] [mi,m2] 

/r C [mia(TOi), m2a(m2)]. (39) 
T/ien (^PrJ admits a radial solution u and 

mi < lr{u) < 7712 a.e in fl. (40) 

For the proof, we refer the reader to 

Generalizing this construction type of the diffusion coefhcient a we obtain 

Proposition 3.4. Let r G [0,d]. Assume that \3'A) holds true and there exist 
an odd integer ni and ni + 1 positive real numbers {?Tii}i=o...ni ; with toq — 
and for all i G {0, . . . , rii — 1} we have mi < m^+i. Moreover 

a{mi) — max a; a(mi+i) = min a Vi G {0, 2, . . . , ni — 3, ni — 1} 

[rai^nii^i] [vai.'mi-^i] 

Ir C Pi [mia(TOi), mj+ia(TOi+i)] 

z— 0.2,...,ni— 3,ni — 1 

(41) 

Then (Pr) admits at least ^^2^ radial solutions {ui}i^^Q2,...,ni-i] such that 
mi < lr(ui) < TTii+i Vi G {0, 2, . . . , rii — 3, rii — 1}. 
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Proof. The proof here is by induction. Indeed we set 

'Pni — { If condition (HT]) is satisfied then(Pr) admits at least ^ — solutions.} 

By using lemnia [?31 with mi — and m2 = mi, it is easy to prove for ni = 1 
that Vni is true. For ni > 1, This procedure can be repeated to prove that if 
'Pni-2 holds true then Vni holds too. □ 



Example 1. Let us see a function a satisfying proposition \S.4\ For this, we 
consider the case rii = 3 and r £ (0, d] . Considering iSS^) and the strong 
maximum principle we get minir > 0. Taking 

maxir a(0) 

with a(0) > and also a decreasing on [0,mi] then we prove lemm diK^ condi- 
tions. 

By repeating this process with m2 > mi and setting 

minIr maxIr 
a[m2) , ma := 2 ■ 



TO2 ' a(m2) 

with a(m3) :— "^"''^ and also a is decreasing on [m2,m3]. This shows the 
existence of such a. 




mj m3 

Figure 1: The case ni = 3 

In the representation of a we have deliberately left, on solid line parts of 
the curve satisfying the conditions of proposition \3^\ and dotted line one without 
constraints. This situation are explain in the figurd\^ 
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Remark 3. As previously announced, the proposition \S.4\ generalizes a certain 
point of view Theorerr A3.ll However it does not accurately determine the exact 
number of solutions of (Pr) and the bifurcation points of branch of solutions. 
We have shown in way to solve this problem by using the linearized problem, 
the principle of comparisons obtained in and the Krein-Rutman theorem. 

3.3 Stable solutions of (P,.) 

Definition 3.5. Given a domain Q. C M", a solution Ur e Hq{VI) of (Pr) is 
stable if: 

V(/) e Hl{n) GuM) ■■= I a{lr{Ur))\\/(l)f - I a{lr{Ur))lr{<j3)VUry<j) > 0. 



(42) 

Definition 3.6. Given u : [0,(i] — )■ H^{Q.), the graph of u is called a (global) 
branch of solutions if 

(i) u(,c{[o,d],Hl{n)), 

(a) u{r) is solution to (Pr) for all r in [0,d]. 

u is called a local branch if it's defined only on a subinterval of[Q,d] with positive 
measure. 

Before concluding this section, we will focus into the case a nonincreasing to 
prove the stability of the global branch of solutions. Assume for all r £ [0, d], 
Ur is a solution to {Pr) and 

< lr{ur){x) < fid for a.e x G J7. (43) 

Assume that there exists a solution /id to (pijl such that 

o-ifJ-d) = min a and a(0) = max a. (44) 

We prove in [3] 

Theorem 3.7. Assume i3^) . i44\ l o-^id holds. Assume in addition 

that a e W^'°° {R) and for some positive constant e, it holds that 



Cl|5|2|/|2|a'U,[-e,Md+e]^^^^ < 1, (45) 



where Ci is a constant dependent to O. Then 
(i) For all r in [0,d], (Pr) possesses a unique radial solution Ur in [uo,Ud\; 
(a) {{r,Ur) : r G [0,(i]} is a branch of solutions without bifurcation point; 
(Hi) it is only global branch of solutions; 

(iv) if in addition, a is nonincreasing on [0,/irf] then r i-^ Ur is nondecreasing. 
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Remark 4. It is very difficult to obtain property (iv) for any a. However when 
a is nonincreasing provide us important information for studying the stability 
of this branch of solutions. 

Corollary 3.8. Let the smallest solution to (Pd)- Assume hS'^) and i f5I]) 
holds true and there exists a solution fXd to iSl]) satisfied ^J^. Assume in 
addition that a G W^'°°{M.), satisfied and for some positive constant e, 
it holds that 

1 

where Ci is a constant dependent to £7. 

Then {(r, it^) : r 6 [0,(i]} is the only global branch of solutions starting to 



C^i|5|2|/|2|a'|oo,[-e,,..+e]7^^ < 1, (46) 



Proof. The fact that {(r, Ur) : r £ [0,(i]} is the only global branch of solutions 
results from theorem l3.7l We will now show that this unique branch of solutions 
is stable and start &t r = dhy u\. For this we consider without loss of generality 
{Pd) admits two solutions u]^ and such that < u^. We denote by /ii and 
Ii2 respectively solutions of (PT|) corresponding to and (see figurej^l). It is 
easy to see that /ii and /i2 satisfied (|33]). 

Assume {{r,Ur) : r G [0,d]} is the only global branch of solutions starting 
to ul. Then we get Ci\g\2\f\2\a'\oo,[~.e..t,2+e]-^(^ < 1- In this case, using 
theorem l3.7l we get (P^) possesses a unique radial solution u,. in [uo,w^] and the 
mapping r u,. is nondecreasing. By continuity of this mapping, we can find a 
To €]0, d[ such that — for a.e x € fl. This means that is a solution 
of (Pro)- This gives us an absurdity and concludes the proof. □ 

We are now able to prove: 

Proposition 3.9. Under assumptions and notation of corollaru [3.8l the global 
branch of solutions described in theorerr AS. 71 is composed by stable solutions. 

Proof. For all r G [0,d], let Ur be a solution belonging to the global branch of 
solutions described in theorem l3.7l By using the linearized problem of (Pr), we 
get V0 G H^{n) 



a{lr{Ur))\V(j)\^- / a'(lr{Ur))lr{(f>)VUrWcl)> 



n n 



(47) 



Taking into account that |Vur|2 < ^^^^ mia\\^\u )) ^^^^""^ C{Q) designed the 



Poincare Sobolev constant. We obtain 



inia{lr{Ur))\\7(j)\l ~ C \g\2\a'\oo,[-e,f,i+e]NUr\2\'^(l)\l- 

\f\2 

a{lr{Ur)) 

a[lr{Ur))\V(t)\^~ j a'{lr{Ur))lr{(j))VUrV(f> > 

n 

\\7(l)\l(^Ma{lr{Ur)) - Ci \g\2\a\oo,[-e,f,i+e 



I/I2 



ini a(lr(ur)) 



(48) 
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Moreover by assumptions and we get a(/ii) < ini a{lr(ur)). 
Thus (|46)) becomes 

Ci|5|2|/|2|a'U,[-e,Md+£l—r -7— -72 < 1- 

mi allriUr}] 

n 

We deduces 

J a{lriUr))\V(t)\'^ - J a' {lr{Ur))lr{4>)VUrV (j) > 0. 

n n 
This concluded the proof. 



a(0) 

a(\mu_l) 



Figure 2: case of 2 solutions 





"1 



r=0 



Figure 3: Branch of solutions 
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4 Parabolic problem 

4.1 L°° estimate 

In what follows we obtain L°° estimate of the solution ([T]) from L'' estimate. 
The method we use is based on iterations Moser type, for more details on the 
method see |14j. 
We get 

Theorem 4.1. Let n > 3 and u a classical solution of (QP defined on [0,T). 
Assume that p > 1 and q > 1 such that ^ + ^ = 1- Suppose further that 
Ug = supt<y \u{t)\q <oo, f e L°°{0, oo, L'^in)). Ifp < ^ then Uoc < oo. 

To prove this theorem we need the following proposition: 



Lemma 4.2. Consider u a classical solution of (QP on [0,r), r > 1 and p > 

1 such that - + - — 
p 1 

supt<^T\'u.{t)\r} and let 



1 such that p + I — 1 with p < We take Ur = max{l, |wo|oo, Ur 



pin + 2) 

cr(r) - 



2[r{2p ~ pn + n) + up] 
Then there exists a constant C'2 — C'2(S^,m) such that 

< [C2 ||/|U=o(o,«,,i,(o))]"(''V-W[7,. 

Proof. Multiplying ([T|) by u^'''^^ and then using the Holder inequality yields 

^jj u^'dx + J \V{unfdx < \fU\u^'-X. (51) 

n n 

As 

k"-'lp=KI,'&, (52) 
by taking w = in ([5T|) and ([5^ . we get easily 

^|l^l2 + rr>^^^\^Ml < l/U^ISp, (53) 



with a = ■ Let B such that 



ap I 

with 2* = We claim that /3 G (0, 1). 
In fact 

2nr - (n - 2)(2r - \)p 
^ ~ (n + 2)(2r- \)p ' 

Sincep < 57^1 -^hi then /3 > 0. As well as (n+2)(2r-l)p > 2m'-(n-2)(2r-l)p 
implies that /3 < 1 this prove that /3 e (0, 1). 

Using an interpolation inequality (see[14j) in ((53|) and (|54)) . we get 



2r"di'^'^ + "'^^'^^'^-'-^'«('^'i'"''^* ) ■ ^^^^ 
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Applying Sobolev injections in (j55p . we have 



1 1 i9 2?' — 1 ,„ ,9 



and also 

Id,,, 2r — 1 ,„ 



^ m J 



(56) 



2r /' 



(57) 

Since /3 G (0, 1) and f e (0, 1) it is clear that e (0, 1). Applying Young's 

A^e obtain 

a(l-/?) 



inequality to ^ with + 1 - ^^^3^ = 1. We obtain 



1 I i9 2r — 1 ,^ ,0 



I/IK- 



with (5=1 — 



Joining the fact that "^"^^ G (0, 1) and 5 < 1 to we deduce 



1 d 



3r - 2 



2r 



We set 



2r2 



a(l - /?) . ^ /3a 

2ra(r) - 1 = and 2p(r) = ^, 



2^)1^^ 
(58) 



(59) 



((5^ becomes 
1 d 



3r-2 



|Vz«|^<|/||(- 

' TO 



2ro-(r)-l 



(60) 



This gives us taking into account that ^ > 1 



2r\ 



2r(7(r) 



+ m|Vu;|^ < |/|| (^^J to Ct . 

By a calculation we can verify that 

2nr - (n - 2)(2r - l)p 



(61) 



p{r) 



2r{p{n + 2) + n) - 2n(2r - l)p' 



and also that p{r) G (0, 1). 

Using the Poincare Sobolev inequality and that p{r) < 1 in (pTj) . yields 



i /2r\ 



2r(T(r) 

where Ci (SI) designed the Poincare Sobolev constant. Noticing that 



(62) 



dt 



Ci(l)) 



2r\ 



2ra-{r) 



^«|^<|/||(-j \w\imC^. 

(63) 
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and integrating ([55)1 on [0,i) we get 

/ c) \ 2r(T{r) 

Ht)\l < \wml + \\f\\U^o,^^,.^n)) [-^) mC^\w\l 



(64) 



Since 

\w{Q)\j = / wiOfdx = / uiO^dx < \n\\u{Q)\^^ < |^)|C/2^ (65) 



and gives us 

2r(T(r) 



1 /2r\ 
C/|; < + ||/|ll.(o,oo,L,(o)) [-) mCmr. (66) 

Whereas | > 1, 2ra{r) > and cr(r) = ^ it follows that 

C/2.<C2^^'-^||/|lli;o,oo,L.(^^))'^'*''^^- (67) 
with C2 = 6*2(^1, to). This completes the proof of Lemma. □ 
We have also 

Lemma 4.3. Let r > 1, n > 3, p < and a{r) = 2[r{2j'-pn+n)+np\ ^'^^'^ 
get 

cr(2'=r) < 6l'^cr(r) Vfc G N, 

TOf/i 6* e (0, 1). 

Proof. By asking ci = p("+2) ^ _ _ 7^) and C3 = np yields a{r) — 

with ci, C2, C3 € R!j^ .By taking 6* = 1 — ^^^^jj^^^ the proof of this lemma is 
deduced by reasoning by induction. □ 

Returning now to the proof of the theorem, by lemma H?^ we get 

Proof. 

C/2. < [C2 ||/|Uoe(o,oo,L.(n))]"'''V^(^)[/,.. 

By iterating this equation by taking r = h,r = 2h, r = 2^h, etc, we obtain 
with 

Ai := a{h) + a{2h) + a{2^h) + .. + a{2''-^h) + a{2^r) 

et 

A2 a{2h) + 2a{2^h) + Za{2^K) + ... + (fc - \)a[2^-^h) + ka{2^r). 

To complete the proof we just need to show that Ai,A2 < +00. Indeed by 
lemma lT^ 

k 00 , ,\ 

Ai < y a^G(h) < y a^'aih) = \ < 00. 
^-^ ^-^ (1 — a) 

Noting also that 

a{2^h) < d''-^a{2h) V/c G N*, 
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it follows 

tj.=i fj.=i ^ ' 

This completes the proof of the theorem. □ 



4.2 Uniform estimate in time 

We prove in what follows an estimate of u in i°°(R"'", iJg (il)). We get 

Theorem 4.4. Assume that f e L'^{^), g e H^{i}), uq G Hq{V,) and a S 
W'^'°°{R) withinia > 0. Then a solution u of is such thatue L°°(M.+ , H^in)). 

Taking a spectral basis related to the Laplace operator in the Galerkin ap- 
proximation (see |16j) we find that —Am can be regarded as test function in 
1/2(0, T, 1/2 (Q)) for all T > 0. By multiplying ^ by -Au{t) and integrating 
over n 

{ut, -Au) + {-div{a{lr{u))Vu), -Au) = (/, -Am), (68) 

and also 

l^Mlv + {-ailr{u))Au, ~Au) + {-a'{lr{u))Vlriu).Vu, -Au) = if, -Au). 

(69) 

Here (., .) is the usual scalar product on L'^{n). Taking to account 

|V/,(m)|2 <i^||.g||Hi(o)|VM|2, (70) 
where K is a constant depending of f2. It comes 

\{-a'{lr.iu))\7l^{u).\7u,-Au)\ < K \\g\\Hiin)\a'\ooMl.\Au\2 (71) 
Now from (|7ip and ([55)1 we have 

i|||Mf^ + m|AM|2-X||5||^,(,,)|a'|oo||ufy|AM|2< I/I2IAMI2. (72) 

By using Young's inequality ab < -^a? + yfe^, we get 

1 , .,0 1 
at mm 
In order to apply the uniform Gronwall lemma to (j73p we start with a small 
estimate. Remember that 



WWl < -\.f\l + -{K \\9\\m(n))^a'\\U\nt. (73) 



^^\u\l+m\\u\\l<^\f\l (74) 

where A is the principal eigenvalue of the Laplacian operator with Dirichlet 
boundary conditions. 

By integrating on [i,to) we get 

t+to t + to 

\uit + to)\l+m f \\ufyds< f J-\f\lds+\u{t)\l (75) 



m 
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and also 

t+to 



7 " " " A m 

t 



uryds<-^\f\lds + -\uit)\l (76) 



Let po > such that |u(i)|2 < Po- -1^7 setting 

m m TO^ m 

we obtain by using uniform Gronwah lemma to (|73p 

||M(i + <o)lk < (^ + a2)exp(ai) Vt > 0, io > 0. (77) 
to 

Hence u G L°°{to, +oo, ffp (51)). By using ([75]) and the classical Gronwall lemma 
it is easy to see that u G L°°(0, -f^o (^^))- This completes the proof of the 
theorem. 

Remark 5. This theorem show us the existence of absorbing set in HQ{fl). By 
considering S(t) the semigroup associated to the equation (QP defined by 

S{t) : L^{n) L^{n) 

with u{t) a solution of (Qp. As a result of the theorerr ^4-4\ '"^'^ compact 
embedding of HQ{fl) into L^{fl) we deduce that the semigroup S{t) possesses 
a global attractor. Indeed it is easy to show the existence of absorbing set in 
L'^{fl), the main difficulty here is to show that S{t) is such that 

MB C L^{VL) bounded, 3t^ = to{B) 

such that (JS'(t)i? is relatively compact in i^(f2). C^^-^ 

t>to 

This property known as S{t) is uniformly compact for t large can be proved by 
using theorem^J^ and the compact embedding of Hq{^) into L'^{Q,). 



4.3 Asymptotic behaviour 

In this part we are interested in asymptotic behaviour of a weak solution of ([1]) . 
Our main interest here is the radial solutions. By radial solutions we means 
u(t, \x\) — u{t,x). As in the stationary case $7 is a open ball of M". Remember 
that 

Ll{n) = {v e L^ifl) 35 G L2(]0,d/2[) such that v{x) = v{\\x\\)}. 

In order to not make confusion between uq the solution to (Pq) and the initial 
value of ((Ij, we will take u'^ the initial value of 

Theorem 4.5. Assume that /, 5 G L'^{Q,), a is a continuous function and the 
assumption checked then (QP admits a radial solution. 

Proof. Let w G L^(0, t, L^(fi)) it is clear that lr{w) is radial and also a{lr{w)). 
Thus by Fr maps L^{0,t, L^{fl)) into itself. The proof now follows by using 
arguments similar to those used in theorem l2.1l □ 
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Assume now 

/,5>0 in n (79) 

and 

uq<u° < Ud, (80) 

with u'^ the initial value to ([1]) and uq and Ud respectively the solution of (Pq) 
and (Pd). 

We can now give a stability result assuming that ([T]) admits a unique solution. 

Theorem 4.6. Assume i 79{} and f,g e L^(51). Let u, Ud and uq respectively 
the solution of U^, (Pd) and (Pd)- If 

uo <u° < Ud, 

then 

Uq < U < Ud Vi. 

Proof. Let 

S^{t\ l{u{s))e[QMud)] Vs<i}- (81) 
It is easy to prove that S contains (see [80)) . By setting 

t* = sup{t \t E S}. (82) 
By continuity of the mapping t — > ld{u{t)). we get 

ld{u{t*))E[QMud)]- (83) 
By using ^ and {Pd) we get in l?(0,t*) 

J^{ud-u, (/.)+ j a{ld{u))V{ud~u)V(j) = - y {a{ld{ud))-a{ld{u)))VudV(l) ^cj) e hI{9). 
u u 

(84) 

Choising cf) ~ {ud — u) , becomes 

^^I(urf-ti)"l2+ J a{ld{ud))\S7{ud-uy\^ = j {a{ld{u))-a{ld{ud)))VudV{ud-u)~ . 
n n 

. . 

Since a is nonincreasing {a{ld{u)) — a(ld{ud)) > < t*) hence proposition l3.2l 
yields 

J {a{ld{u)) - a{ld{ud)))\7udV{ud - uY < 0. (86) 

Thus 

~\iud - u)-\l + a{ld{ud))\Viud ~ u)-\l < (87) 
Applying Poincarre Sobolev inequality we get 

^^l("d - u)-\l + C2\iud - u)-\l < 0, (88) 

this prove 

l{e^^^^\iud-ur\l}<0. (89) 
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Moreover {ud — u) (0) = {ud — u") = it follows that Ud> u V< G [0, t*]. In 
the same way we can also prove uq < u Vi G [Oji"^]. It follows 

uo<u<Ud Vt e [0, t*] (90) 

To finish we just need to prove that t* is very large, this is typically the case. 
Indeed if < oo then 

l{u{t*)) = or Idiud). (91) 
From dini) and (HO]) we deduce 

u{f*) = Uq or u{t*) = Ud- (92) 

Due to the uniqueness of ((IJ, we deduce that t — oo. This shows that 

Uq < U < Ud Vi, 

and achieve the proof. □ 

Remark 6. The fact that \u{t)\2 is not a Lyapunov function that is to say de- 
creases in time, makes very complex the study of certain asymptotic properties 
of our problem. Indeed under our study it is tempting to show that for r fixed 
r e]0,d[ 

u{t)^ul in L^{n), 

where u is the solution of {Ip and u], the solution belonging to the stable global 
branch described previously. A numerical study would be a great contribution to 
try to carry out some of our theoretical intuitions. 
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